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Abstract 


We present a class of models in the framework of gauge mediation supersymmetry 
breaking where the standard model is supplemented by additional U(l) symmetry 
which acts only on the third generation fermions. The messenger fields carry non¬ 
trivial U(l) charge and are vector-like particles under this symmetry. This leads 
to additional contribution to the soft supersymmetry breaking mass terms for the 
third generation squarks and sleptons. In this framework we show that the muon 
g — 2 anomaly, the observed 125 GeV Higgs boson mass and the detected relic dark 
matter abundance (gravitino in our case) can be simultaneously accommodated. The 
resolution of the muon g — 2 anomaly, in particular, yields the result that the first 
two generation squark masses, as well the gluino mass, should be < 2.5 TeV, which 
will be tested at LHC14. 
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1 Introduction 


The discovery of the Higgs boson with a mass ~ 125 GeV [1, 2] has important conse¬ 
quences in general for low scale supersymmetry (SUSY). In the minimal supersymmetric 
standard model (MSSM) it requires either a large, (9(few — 10) TeV, stop quark mass, or 
alternatively a relatively large soft supersymmetry breaking (SSB) trilinear H^-term, along 
with a stop quark mass of around a TeV [3]. The constraints related to the Higgs boson 
are particularly stringent regarding the sparticle spectrum in the gauge mediated SUSY 
breaking (GMSB) scenario [4], In both the minimal (mGMSB) [5] and more general [6] 
GMSB scenarios, the trilinear SSB H—terms are relatively small at the messenger scale 
and hence, accommodating the light GP-even Higgs boson mass of 125 GeV requires a 
stop mass in the multi-TeV range [4], This, in turn, pushes the remaining sparticle mass 
spectrum including sleptons to values above the TeV scale or so [4] . Imposing t — b — r (or 
b — t) Yukawa coupling unification condition at Mqut in the GMSB scenario makes the 
sparticle spectrum even heavier than this. It was shown in Ref [7] that this latter scenario 
predicts that all colored sparticle masses are above 3 TeV, which would be hard to test at 
LHG 14 but may be accessible at some future colliders such as HE-LHG 33 TeV or a 100 
TeV collider. 

The Standard Model (SM) prediction for the anomalous magnetic moment of the muon, 

= {g — 2)^/2 (muon g — 2) [8], has a (2 — 3) a discrepancy with the experimental results 
[9]: 

Aa^ = a^(exp) — a^(SM) = (28.6 ± 8.0) x 10“^° . (1) 

If supersymmetry is to offer a solution to this discrepancy, the smuon and gaugino (bino or 
wino) should be relatively light, 0(100) GeV. Thus, it is hard to simultaneously explain the 
observed Higgs boson mass and resolve the muon g — 2 anomaly with mGMSB (mGMSB). 
There have been several recent attempts to reconcile this (presumed) tension between muon 
g — 2 and 125 GeV Higgs mass in the mGMSB framework. For instance, adding low mass 
vector-like particles [10] can provide signihcant contributions to the Higgs boson mass [11]. 
Another suggestion is to introduce matter-messenger mixing [12], but this can potentially 
reintroduce the SUSY flavor problem in mGMSB, which is undesirable. Ref. [12] assumes 
a perfect alignment of appropriate parameters in order to avoid unwanted flavor violating 
processes. Alternatively, a specific flavor symmetry is employed to overcome this challenge 
[13]. It was shown in Ref. [14] that by enlarging the messenger sector with a weak triplet 
and a color octet, one can resolve the muon g — 2 anomaly and obtain a 125 GeV Higgs 
boson in the spectrum. In order to satisfy the cosmological and accelerators constraints in 
this model, R-parity violating interactions are introduced. 

Recently, Ref. [15] proposed a class of supersymmetric models in the framework of 
gravity mediated supersymmetry breaking [16], in which symmetry considerations alone 
dictate the form of the soft SUSY breaking (SSB) mass terms for the sfermions [15]. It is 
called flavor symmetry based MSSM (sMSSM), and in this framework the hrst two family 
sfermion masses have degenerate masses at Mqut, while the SSB mass term for third family 
sfermions is different. It was shown in [17] that constraints from flavor changing neutral 
current (FGNG) processes, for the case when third generation sfermion masses are split 
from masses of the first and second generations, are very mild and easily satisfied. This 
approach therefore allows for significantly lighter first two family sfermions, while keeping 
the third generation sfermions relatively heavy. The hrst two families and gauginos (bino 
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and wino) play an important role in the resolntion of the mnon g — 2 problem. The anthors 
in Ref. [17, 18] assnme MSSM Higgs fields {Hu and Ha) with a common SSB mass either 
with the first two generation sfermions or with the third generation sqnarks and sleptons. 
In these approaches it was shown [18] that one can easily accommodate the resolntion 
of mnon g — 2 problem with the 125 GeV Higgs boson mass, bnt it is difhcnlt to obtain 
nentralino dark matter with the correct relic abnndance. Ref. [19] proposed to consider 
SSB mass terms for Hu and Hd independent from both the hrst two generation and third 
generation sfermion SSB mass terms. As a conseqnence, it is possible to simnltaneonsly 
have 125 GeV Higgs boson, solve the mnon g — 2 anomaly and accommodate the correct 
relic abnndance using nentralino- slepton coannihilation. Moreover, this framework allows 
one to implement t-b-r Yukawa unification [20]. 

In this paper we propose two models which allow us to realize the sMSSM structure 
in the GMSB framework, with relatively light and nearly degenerate squark and slepton 
masses of the first two families. The third generation sfermions are heavy with masses 
~ 10 TeV. One simple way to realize this scenario within GMSB is to introduce a flavor 
dependent U{1) gauge symmetry, and we present two distinct scenarios here. We hrst 
consider f/(l)i, under which all the left-handed fermions from the third generation have 
charge i+'g), while the third family right-handed fermions (including right-handed tan 
neutrino) have charge {-rj)- We also assume that the messenger helds (for instance (5 -|- 
5)) transform non-trivially under f/(l)i. Thus, the messenger helds interact with the 
third generation sfermions via U{l)i gauge interaction and provide additional contributions 
to their SSB mass terms. This new contribution makes the third generation sfermions 
adequately heavy, while the hrst two family sfermions remain light and degenerate in mass. 
Note that we also increase the contribution to the trilinear SSB breaking terms for third 
generation sfermions through the U{l)i gaugino loop, but this additinoal contribution turns 
out to be negligible. The model allows mixing between the hrst two generation fermions 
with the third generation via a non-renormalizable coupling. However, the desired mixing 
between the generations which yields the correct Gabibbo-Kobayashi-Maskawa (GKM) 
mixing matrix is realized after spontaneous breaking of f/(l)i. 

In Model H a f/(l )2 symmetry acts only on the third generation fermions. In this case 
the left-handed top and bottom helds have a U{1)2 charge of + 7 ], while the left-handed 
tan lepton charge is -Sg. The right-handed top and bottom helds have charge —g, and 
the right-handed tan and tan neutrino have charge +3g. Glearly, in Model H the U{1)2 
charge assignment resembles the well-known U{1)b-l assignments, except that here U{ 1)2 
is a havor symmetry which acts only on the third generation fermions. 

This paper is organized as follows: We describe Models I and H in Section 2 and 
in Section 3 we summarize the scanning procedure and the experimental constraints we 
employ. In Section 4 we present our results focusing on the resolution of muon g — 2 and 
accommodating the 125 GeV Higgs boson mass and relic dark matter abundance. We also 
provide two tables of benchmark points in this section, which exemplify our hndings. Our 
conclusions are presented in Section 5. 


2 Essential Features of the Model 

Supersymmetry breaking in a typical GMSB scenario takes place in a hidden sector, and 
this effect is communicated to the visible sector via messenger fields. The latter interact 
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with the visible sector via the SM gauge interactions, and induce the SSB terms in the 
MSSM through loops. In order to preserve perturpative gauge coupling unification, the 
minimal GMSB scenario can include (5 + 5) {N^ = 1,...,5) SU{5) multiplets or one 
(10 +10) pair, or 10 + 10 + 5 + 5, or one pair of 15 +15. For simplicity, we only consider the 
case with (5 + 5) vectorlike multiplets. Also it is known [4] that the sparticle spectrum 
does not change drastically for > 1 compared to the A^s = 1 case. The (5 + 5) includes 
SU{2 )l doublets (L + Z), and SU{3)c triplets {qc + Qc)- Note that the (5 + 5) fields carry 
additional U{l)i (or U{1)2) charges {rj + (—//)). In order to incorporate SUSY breaking 
in the messenger sector, the fields in (5 + 5) multiplets are coupled, say, with the hidden 
sector gauge singlet chiral field S: 

WDyiSii + y 2 Sqq, (2) 

where W denotes the appropriate superpotential. Assuming non-zero vacuum expectation 
values (VEVs) for the scalar and F components of S, namely S = {S) + 9‘^{F), the mass 
spectrum of the messenger fields is as follows: 

TJif = M. (3) 

Here mb and m/ denote the masses of the bosonic and fermionic components of the appro¬ 
priate messenger superfield, M = y{S) and A = {F)/{S). The dimensionless parameter 
A/M determines the mass splitting between the scalars and fermions in the messenger 
multiplets. This breaking is transmitted to the MSSM particles via loop corrections. 

At the messenger scale we have the gauge symmetry SU{3)c x SU{ 2 ) lxU (l)y x f/(l)i 
(or 17(1)2), where the 17(1) charge assignments are listed in Table 1. The gaugino masses 
are generated at 1-loop level, and assuming {F) -C {S)^, are given by 

rv • 

M, = N,A-^, (4) 

47r 

where i = 1, 2, 3, 4 stand for the SU{3)c, SU{ 2 )l, U{1)y and 17(l)i (or 17(1)2) sectors, 
respectively. The MSSM scalar masses are induced at two loop level: 

i=l 

where Ci = 4/3, C 2 = 3/4, C 3 = (3/5)(Y/2)^ and = 7 ]“^. Here Y and rj denote the 
hypercharge and 17(l)i (or 17(1)2) charges respectively. 
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Table 1: 17(1) charge assignments for the third generation fermions. Here we use the 

standard notation for the SM fermions. 17(l)y stands for the hypercharge and 17(1) 1 (and 
17(1)2) are additional gauge symmetries that we introduce for Model I (and Model H). 

As seen from Eq. (5), the extra contribution to the sfermion SSB masses from 17(1) 1 
(or 17(1)2) gauge mediation is proportional to (74 = and 0:4. A suitable choice of these 
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parameters can make the contribution from U{l)i (or [ 7 ( 1 ) 2 ) gauge mediation independent 
from what we have in mGMSB. Of course, this contribution cannot be completely arbi¬ 
trary, but we assume it to be ten times or so larger compared to the contribution from 
SU{3)c X SU{2)l X [7(l)y. Note that the [7(l)i (or [ 7 ( 1 ) 2 ) contribution, denoted as D, in 
principle, can be quite independent from the SU{3)c x SU{2 )l x U{1)y sector if we intro¬ 
duce messenger helds with only [7(l)i (or [[( 1 ) 2 ) charges, and which are singlets under the 
SM gauge symmetry. In this case we do not need to worry about = 77 ^ and values 
at all. We can just use proper values of A' in Eqs. (5) and (6). Since the origin of new 
contributions to the third generation sfermions can not be tested at the experiment we do 
not specify the model and added just a new parameter D at the messenger scale which 
measures the new contribution to the third family. 

The summary of model building part of the scenario. New contribution to the third 
generation sfermions are obtained by extending the gauge symmetry by an additional U(l) 
symmetry compared to the canonical MSSM gauge mediation scenario. Note that based on 
the charge assignment presented in Table 1 the Model I has a chiral anomaly under the new 
[7(l)i symmetry. This anomaly can be canceled by introducing set of new particles. For 
simplicity, we assume that the [7(1) 1 symmetry is broken not far below the messenger scale 
and these extra particles introduced for just cancellation of the chiral anomaly, and hence 
they obtain masses around the [7(1) 1 symmetry breaking scale. So below the messenger 
scale RGE evaluation is exactly the same as what we have in MSSM. Existence of a new 
U(l) symmetry in our model changes only the boundary conditions on the SSB mass terms 
for the third generation. Since there are a number of ways to cancel [7(l)i anomalies and 
the new particles are superheavy we do not specify them here. 

On the other hand the situation is much simpler for Model II case. The [7(1)2 charge 
assignments resemble U{1)b-l charges (see Table 1) except that the former is a flavor 
symmetry which acts only on the third generation fermions. In this context. Model II is 
very minimal in terms of additional particles in comparison to the canonical MSSM gauge 
mediation scenario. We just need to introduce the held which can break spontaneously 
£/(1)b-i symmetry. Again here for simplicity, we assume that [7(1) b-l symmetry is broken 
just below the messenger scale; therefore, the RGE is the same as which we have in the 
MSSM case. 

The A-terms in GMSB models vanish at the messenger scale (except when the MSSM 
and messenger helds mix [22], which we do not consider in this study). They are generated 
from the RGE running below the messenger scale. The sparticle spectrum in our model(s) 
is therefore completely specihed by the following parameters dehned at the messenger scale: 


ATmess) A, tan/3, S7^?7,(/i), A^5, Cgrav) D. 


( 6 ) 


Here M^ess = M and A are the messenger and SSB mass scales dehned above , and tan/3 
is the ratio of the VEVs of the two MSSM Higgs doublets. The magnitude of /i, but not 
its sign, is determined by the radiative electroweak breaking (REWSB) condition. The 
parameter Cgrav(> 1) ahects the mass of the gravitino and we set it equal to unity from 
now on. For simplicity, we consider the case = 1. Ghanging the value of does not 
signihcantly alter the sparticle spectrum [4]. 

In summary, the MSSM sfermion masses have the following expression: 


= ("^0i,2)gmsb, rnl^ = (m|3)GMSB + 


(7) 
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where 0 i ,2 and 03 denote the sfermions of the first two families and the third family 
respectively, and stands for the sparticle charges under 17(1)i (and f/(l) 2 ). We can 
normalize the sparticle charges by setting rj = 1. 

3 Scanning Procedure and Experimental Constraints 

We employ the ISAJET 7.84 package [27] to perform random scans over the fundamental 
parameter space. In this package, the weak scale values of gauge and third generation 
Yukawa couplings are evolved to Mmess via the MSSM RGEs in the DR regularization 
scheme.The various boundary conditions are imposed at Mmess and all the SSB parameters, 
along with the gauge and Yukawa couplings, are evolved back to the weak scale Mz. 
The threshold corrections [28] are added to the Yukawa couplings at the common scale 
AfsusY = where and are the soft masses of the third generation left 

and right-handed stop quarks respectively. The entire parameter set is iteratively run 
between Mz and M^ess, using the full 2-loop RGEs, until a stable solution is obtained. 
To better account for leading-log corrections, one-loop step beta functions are adopted for 
the gauge and Yukawa couplings, and the SSB parameters m* are extracted from RGEs at 
multiple scales: at the scale of its own mass, m, = for unmixed sparticles, and the 

common scale Mgugy for the mixed ones. The RGE-improved one-loop effective potential 
is minimized at Mgusv, which effectively accounts for the leading 2-loop corrections. Full 
1 -loop radiative corrections [28] are incorporated for all sparticle masses. 

We have performed random scans over the model parameters given in Eq.(6) in the 
following range: 


10^ GeV < 

A 

< 10^ GeV, 

A GeV < 

Afmess 

< 10^® GeV, 

0 < 

D 

< 30 TeV 

2 < 

tan (3 

< 60, 

N 5 = 1, 

/i > 0 

rrit = 173.3 GeV. 


W can be varied from 1 to 4, but the sparticle spectrum does not change significantly [4], 
thus we set W = 1 for entire scan of the parameter space. Regarding the MSSM parameter 
p, its magnitude but not its sign is determined by the radiative electroweak symmetry 
breaking (REWSB). In our model, we set sgn{fi) = 1. Finally, we employ the current 
central value for the top quark mass, rrit = 173.3 GeV. Our results are not too sensitive to 
one or two sigma variation of m* [29] . 

In scanning the parameter space, we employ the Metropolis-Bastings algorithm as de¬ 
scribed in Ref. [30]. The data points collected all satisfy the requirement of REWSB. 
After collecting the data, we impose the mass bounds on the particles [31] and use the 
IsaTools package [32] to implement the various phenomenological constraints. We succes¬ 
sively apply mass bounds including the Higgs [1, 2] and gluino masses [33]. In addition, 
we apply the mass bounds on selectrom/smuon as > 245 GeV if the smuon is NLSP 
[34] and decays inside the LHG detector. But if the selectrom/smuon are long live enough 
in order to decay outside of detector then > 400 GeV [42]. Beside the mass bounds, 
we also apply the constraints from the rare decay processes Bs —>■ [36], b ^ sj [37] 
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and Bu —)■ tUt [38]. The constraints are snmmarized below: 


123 GeV <mh< 127 GeV 
rrig (and m^) > 1.5 TeV 

> 245 (400) GeV (if smnon is NLSP) 


0.8 X 10"^ < BR{Bs < 6.2 x 10"® (2a) 

2.99 X 10“^ < BR(6 ^ 57 ) < 3.87 x 10“^ (2a) 


0.15 < 


BR(i?u —)■ TUr 


IMSSM 


BR(R„ -)■ Tl/-r)sM 


< 2.41 (3a) . 


( 9 ) 


4 Results 


4.1 Model I 

In this section we present the resnlts of the scan over the parameter space listed in Eq. (8) 
for Model I. As previonsly mentioned, the characteristic fntnre of Model I is that all 
third generation sfermions receive additional nniversal contribntion to their SSB mass 
compared to the hrst generation. Fignre 1 represents plots in M^ess — A, Mmsss — h. A — D, 
and Mjness ~ D planes. All points are consistent with REWSB. Green points satisfy the 
mass bonnds and constraints from rare B-meson decays. Yellow points form a snbset of 
green ones, and they represent valnes of Aa^ that wonld bring theory and experiment in 
agreement to within la. We see from the Mmess — A panel that the resolntion of mnon g — 2 
problem allows only a relatively narrow range for A, 5.2 < log (A/GeV) <5.6, while it is 
possible to hnd solntions for 5 < log(Mmess/GeV) < 14. Similarly the /i— term can lie in a 
wide range ~ 2 —10 TeV consistently with all the experimental constraints and compatible 
with the mnon g — 2. The contribntion from the flavor U{l)i symmetry, qnantihed by D, 
is bonnded from the mnon g — 2 constraint (see yellow points). As seen from the Mmess — D 
panel, in order to stay within a la range of mnon g — 2 and have a light GP even Higgs 
boson mass aronnd 125 GeV, the additional contribntion from flavor f/(l)i satishes 5 TeV 
< D < 20 TeV. 

In the minimal GMSB scenario, at the messenger scale the bino and right-handed 
slepton masses are related to each other, namely = |m| [41]. Even thongh all right- 
handed sleptons are degenerate at the messenger scale, the lightest stan can be the lightest 
slepton at low scale dne to the Yukawa coupling contribution to the RGE evolution. Also, 
there is the tri-linear SSB coupling contribution to the stan mass matrix which can make it 
much lighter. Hence, in GMSB we have either a stan or neutralino NLSP depending on the 
initial parameters. However, in our model, since we have additional positive contributions 
to the third generation sfermion masses the stan becomes heavier and so the right-handed 
smnon or selectron can be the NLSP. Note that a right-handed selectron/smuon NLSP 
in gauge mediation models was studied in Ref. [42] without addressing the muon g — 2 
anomaly. Also the sparticle spectrum are very different in comparison since we have a very 
different gauge sector. 

Figure 2 displays the results in the m^o — m^o — mp,^ planes. The color coding is 
the same as in Figure 1 except that the smnon and neutralino mass bounds are not applied 
here. The unit line indicates that it is possible to have neutralino - selectron/smuon 
degeneracy. Note that we have either smnon or neutralino NLSP solutions in the model. 
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log(A/,„oss/GeV) log(il/„.ess/Gc3V) 





//(TeV) 


Figure 1: Plots in M^ess — A, Mmbss — A — D, and Mmess — D planes. All points are 
consistent with REWSB. Green points satisfy the mass bounds and constraints from the 
rare B-meson decays. Yellow points form a subset of green ones, and they represent values 
of Aa^ that would bring theory and experiment in agreement to within la. 


As seen from the m^o — panel, there are plenty of solutions for both neutralino NLSP 
and right-handed smuon NLSP. Even though the neutralino mass bound is not applied in 
this panel, the lightest neutralino mass cannot be lighter than 250 GeV due to the current 
mass bound on gluino {rrig > 1.5 TeV). The reason why the gluino and neutralino masses 
are tied up is due to the following relation among the gauginos at the messenger scale: 
Ml : M 2 : M 3 = q;i : 0:2 : 0 : 3 . Similarly, as mentioned above, the neutralino and the 
right-handed smuon masses are related to each other (M^ = |m| ) as well. However, the 
RGE evolution can lower the right-handed smuon mass to ~ 100 GeV, while the neutralino 
cannot be lighter than 250 GeV. On the other hand, such light smuons can be excluded by 
the search for prompt decay inside the detector if gravitino is the LSP. In this case, the mass 
bound on the right-handed smuon NLSP is > 245 GeV, and it becomes more severe if 
the smuon (and selectron) can decay outside the detector (m^ > 400 GeV) [42]. All these 
bounds are applicable if we assume R-parity conservation. If we assume phenomenologically 
acceptable R-parity violation in the model, then this bound will disappear and one can 
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Figure 2: Plots in m^o m^o planes. The color coding is the same as in Figure 

1 except that the mass bound on smuon and neutralino are not applied here. The unit 
line indicates the regions where the two masses represented on the axes are equal to each 
other. 

have rrifi > 100 GeV. Note that the authors in Ref. [42] presented several models and 
performed the collider phenomenology with smuon/selectron as NLSP. Our model has a 
very different sparticle spectrum compared to the models in Ref. [42] except for the NLSP 
smuon / Solectron. 

The m^o — panel in Figure 2 also indicates that we have a stringent lower bound 
for left-handed smuons/selectrons. As for the right-handed smuon and neutralino, the left- 
handed smuon mass is related to the wino mass as M| ~ f^e messenger scale. 

Since the gluino mass bound affects the wino mass, it also affects the left-handed smuon 
mass. This can be seen from the m^o — rrin, panel in which we only have solutions with 
> 500 GeV. 

In Figure 3 we illustrate our results in the Aa^ — Mmess and — planes. Gontrary 
to previous figures, the gray points are consistent with REWSB and represent solutions 
with a selectron/smuon NLSP. Green points satisfy the mass bounds and constraints from 
the rare B-meson decays. Yellow points form a subset of green ones, and they represent 
values of Aa^ that would bring theory and experiment in agreement to within la. The 
Aa^—Mmess plane shows that we have the desired contribution to the muon g — 2 calculation 
from supersymmetric sparticles for almost any value of Mmess with smuon NLSP. On the 
other hand, this broad range of M^^ss provides a very wide interval for the gravitino mass 
which, in gauge mediation, can be expressed as 


UIq ~ 


T/mess A 
■\/3 Mpi 


( 10 ) 


where Mpi is the Planck mass. As shown in the M^ess ~ A plane of Figure 1, the solutions 
corresponding to the experimental measurement for muon g — 2 within la deviation are 
located in the very narrow interval 1.2 x 10^ GeV < A < 4 x lO^GeV. Having solutions in 
a such narrow interval for A makes the gravitino mass almost linearly dependent on the 
scale Afjijess- 
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Figure 3: Plots in the Aa^ — Mmess and uiq — rrijXj^ planes. In contrary to other hgures the 
gray points are consistent with REWSB with the selectron/smuon as NLSP. Green points 
satisfy the mass bounds and constraints from the rare B-meson decays. Yellow points form 
a subset of green ones, and they represent values of Aa^ that would bring theory and 
experiment to within la. 
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In the second plot in Figure 3 we show the relation between rriQ and with smuon 
being the NLSP. This relationship is very important, as mentioned above, in estimating 
the smuon life-time which decays into a lepton and a gravitino. In general, for the NLSP 
slepton (iji) case the two-body decay width is given by [43]: 


r(4 ^ iG) 



167r(il/messA)2 



( 11 ) 


Using Eq.(ll) and following the discussion in Ref. [42] we show our results in the — 
plane with a solid and a dashed line. The solid line corresponds to 1 mm free length of the 
smuon, while the dashed line to 10 m which is supposedly the length of the detector. If a 
prompt decay is identihed with the 1 mm path in which the smuons fly before decaying to 
a gravitino and a lepton, the region between the solid and dashed lines can be identihed as 
smuon decaying inside the detector. According to the results represented in the ttIq — rrijXn 
plane. Model I predicts a smuon NLSP of mass upto 400 GeV. Hence, the model predicts 
plenty of solutions testable inside the detectors, even if one assumes R-parity conservation 
which excludes the region with < 245 GeV. 

In Figure 4 we show the results in Aa^ — tan/3, Aa^ — Aa^ — and 

Aa^ — rrih planes. The color coding is the same as in Figure 1 except that yellow points 
form an independent subset of gray, and they represents the values of Aa^ which would 
bring theory and experiment into agreement within la. Besides, the Higgs boson mass 
bound is not applied in the Aa^ — mt panel. The Aa^ — tan (3 plane shows the tan {3 
dependence of muon g — 2. Since the contribution from the supersymmetric particles to 
the muon g — 2 depends linearly on tan (3 it is understandable why there are more solutions 
satisfying the muon g — 2 constraint with increasing tan (3. According to the results in this 
plane, it is hard to have substantial contribution to muon — 2 if tan (3 < 30. 

The Attfj, — rurj^/rriej^ and Aa^ — mi^/rriq^ planes highlight the mass splitting necessary 
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Figure 4: Plots in Aa^ —tan/?, Aa^ — Aa^ — mi^/rriq^, and Aa^ — rrih planes. The 
color coding is the same as in Figure 1 except that yellow points form an independent subset 
of gray, and they represent the values of Aa^ which would bring theory and experiment to 
within la. Besides, the Higgs boson mass bound is not applied in the Aa^ — ruh plane. 


between the third and hrst two-family sfermion SSB masses in order to satisfy all current 
phenomenological constraints including muon g—2. As seen from the Aa^—plane, 
the muon g — 2 solution corresponds to the slepton splitting requirement ruf^/mg^ > 8 
at low scale. At the same time for the squarks we have a relatively milder constraint 
as 3 < mi^/rriq^ < 10. Finally, the Aa^ — rrih plane shows plenty of solutions which 
simultaneously accommodate the 125 GeV Higgs boson and muon g — 2 anomaly. 

We summarize the sparticle spectrum in Figure 5 with plots in — tan /?, — rrig, 

mi^ — rriuji, ~ planes. The color coding is the same as in Figure 1. The 

tua — tan (3 plane shows that niA can be as low as 2 TeV for tan (3 > 30, which is required 
also by muon g — 2. It is very interesting to note that the solutions which accommodate the 
muon g — 2 anomaly at la level predict squarks and gluino masses lighter than 3.5 TeV. 
Furthermore, as seen from the — 'n^g plane most of the solutions are associated with 
the gluino masses less than 3 TeV. This makes our model readily testable at LHC Run H. 
The lightest stop is found to be relatively heavy in order to accommodate the 125 GeV 
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Figure 5: Plots in — tan/?, — rrig, mi^ — rriu^, and nif^ — rrig^ planes. The color 
coding is the same as in Figure 1. 


Higgs boson in the model. According to the — rriu^ panel, mi^ > 5.5 TeV for rrtu^ > 2 
TeV. Finally, as seen from the nir^ — plane is bounded in the ~ 200 — 500 GeV 

range as expected, since the supersymmetric contributions to muon g — 2 rely on the light 
smuons. Again, rrif-^ > 5.5 TeV because of the flavor symmetry contribution. 

We present our results for the gravitino mass in Figure 6 with the Aa^ — tUq panel. 
The color coding is the same as in Figure 4. In GMSB scenarios the gravitino is usually 
the LSP, and its mass can vary between 10 TeV and a few eV. As seen from the left panel 
of Figure 6, it is possible to stay within the la band of muon g — 2 for a wide range of 
the gravitino mass, from a few eV to 10 GeV. A light gravitino provides a plausible dark 
matter candidate and it can also manifest itself through missing energy in colliders [39] . In 
standard scenarios, the relic density bound ~ 0.11 [40]) is satisfied with a gravitino 
mass ~ 200 eV [39], which makes the gravitino a hot dark matter candidate. Hot dark 
matter, however, cannot compose more than 15% of the total dark matter density and this, 
in turn, implies that the gravitino mas should be less than 30 eV [39]. The right panel of 
Figure 6 is the same as the left panel, except that it is zoomed to the region which yields 
gravitino masses of the order of a few eV. The predicted gravitino mass compatible with 
muon — 2 is consistent with the hot dark matter gravitino, since resolution of muon g — 2 
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Figure 6: Plot in Aa^ — plane. The color coding is the same as in Figure 4. 


within the sGMSB framework requires > 10 eV as seen from the right panel. In order 
to have a complete dark matter scenario one could invoke axions as cold dark matter in 
this region. 

A gravitino mass > 30 eV requires non-standard scenarios in order to agree with obser¬ 
vations. Such non-standard scenarios include gravitino decoupling and freezing out earlier 
than in the standard scenario, which may be possible in a theory with more degrees of 
freedom than the MSSM [47]. A gravitino of mass > keV is still possible and it can be 
cold enough to constitute all of the dark matter if non-standard scenarios such as early 
decoupling is assumed. 

Finally we present three benchmark points in Table 2 which exemplify our hndings. All 
points are chosen to be consistent with the experimental constraints and muon g — 2. All 
masses and scales are given in GeV units except the gravitino mass which is given in eV. 
Point 1 represents a solution with a the 125 GeV Higgs boson and Aa^ ~ 28.7 x 10“^° with 
the lowest possible messenger scale (M ~ 10®). Point 2 depicts a solution with the least 
amount of mass splitting within the squark and slepton families compatible with the 125 
GeV Higgs and desired muon g — 2 11, mi^/mu^ ^ 4). Points 1 and 2 display 

solutions with neutralino NLSP, while point 3 shows a solution for the smuon NLSP which 
is compatible with the Higgs boson mass, muon g — 2. 

4.2 Model II 

In this section we present our results for a different charge conhguration under the extra 
U (1) symmetry listed as Model H in Table 1. Note that the fundamental parameters remain 
the same, and only the contribution from the extra U{ 1)2 symmetry to the sparticles is 
now different. In Model H, the third family squark masses receive additional contributions, 
as in Model I, parameterized by D^. On the other hand the third family charged slepton 
mass^ receives the contribution of As we will show below, having stau heavier than 
stop has consequences for muon g — 2 contribution and the Higgs boson mass. 

We present our results for Model H hrst in terms of the parameters m^o and which 
are some of the crucial players for muon g — 2 anomaly. The m^o — plot is shown 
in Figure 7 with the color coding the same as in Figure 1. The results displayed in the 
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Point 1 

Point 2 

Point 3 

A 

0.25 X 10® 

0.26 X 10® 

0.34 X 10® 

A/mess 

0.9 X 10® 

0.50 X 10® 

0.20 X lO'^ 

tand 

58 

57 

59 

D 

11580 

9755 

22060 

rrifjms^ 

14 

11 

19 


5 

4 

7 


28.9 X 10-“ 

21 X 10-“ 

23 X 10-“ 

g 

5404 

4424 

10179 

At 

-374 

-371 

-496 

Ah 

-385 

-382 

-509 

Ar 

-36 

-36 

-52 

rrih 

125 

125 

126 

rriH 

2201 

1954 

3504 

rriA 

2187 

1942 

3481 

mH± 

2204 

1957 

3505 

m^o 

Xl,2 

347, 678 

371,723 

483, 936 

m^o 

-^3,4 

5312, 5312 

4347, 4347 

10005, 10005 

m-± 

A1,2 

680, 5245 

724, 4292 

939, 9879 

m-g 

1902 

2012 

2563 


2396, 2279 

2561, 2440 

3060, 2898 

^5.2 

11084, 11332 

9483, 9770 

20940, 21360 

^dL,U 

2397, 2266 

2562, 2427 

3061, 2879 

mt 

01.2 

11275, 11423 

9648, 9770 

21276, 21528 


798 

845 

1031 

mc,^ 

11407 

9642 

21654 


853, 363 

883, 398 

1170, 405 

^n,2 

11067, 11399 

9366, 9637 

20976, 21360 

mg (eV) 

54 

31 

157 


Table 2: Benchmark points for exemplifying our results. All points are chosen to be 
consistent with the experimental constraints and muon g — 2. All masses and scales are 
given in GeV units except for the gravitino mass which is given in eV. Point 1 represents 
a solution with a 125 GeV Higgs boson, Aa^ ^ 28.7 x 10“^° with the lowest possible 
messenger scale (M ~ 10®). Point 2 depicts a solution with the least amount of mass 
splitting within the squark and slepton families compatible with the 125 GeV Higgs and 
desired muon g — 2 {mr^/me^ ^11, mi^/rriu^ ~ 4). Points 1 and 2 display solutions with 
neutralino NLSP, while point 3 shows a solution with smuon NLSP which is compatible 
with the Higgs boson mass, muon g — 2. 
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Figure 7: Plot in the m^o — plane. The color coding is the same as in Figure 1. 

m^o — rrtjXj^ plane exhibit sharp lower bounds on the neutralino and right-handed smuon 
masses. As in Model I, gluino mass bound determines the lower bound on the neutralino 
mass. We have two different bounds for the right-handed smuon mass. If the lightest 
smuon is heavier than the neutralino, the lower mass bound is 100 GeV [31]. When smuon 
is NLSP it decays to gravitino and muon and the lower mass bound 245 GeV. This bound 
can be go up to 400 GeV when gravitino is heavier then keV or so. But this bound can 
be relaxed up to 245 GeV if we assume suitable R-parity violation in the theory. The tail 
shape in the plot corresponding to the region with since the right-handed 

smuon is not NLSP in this region, the ft, —?• fiG decay constraint is not applicable. As we 
can see in this model neutralino can be as light as about 200 GeV or so. There is difference 
compare to the Model I in terms of upper bound for neutralino and smuon. Here the 
neutralino needs to be lighter than 400 GeV and smuon lighter than 320 GeV in order to 
bring the muon g — 2 into agreement within la of the experimental results. It was noticed 
in Ref. [46] that a smuon in this mass interval can be tested at the LHG. On the other 
hand, the bound for this particle in model I was up to 600 GeV. The reduction of the upper 
bound of this particle can be related to the tan /3 bounds for model II which we display 
below. 

We use Figure 8 to illustrate the impact of the extra heavy stau lepton in the spectrum 
and the results are presented in the Aa^ — tan (3, mt — tan m^o — tan /3, and — tan {3 
planes. The color coding is the same as in Figure 1 except for the Aa^ — tan/3 plane 
where the color coding is the same as in Figure 5. The Aa^ — tan f3 plane reveals a sharp 
fall in muon g — 2 values (green points) for tan f3 ~ 40. One would have expected the 
supersymmetric contributions to muon — 2 to linearly increase with tan /3 [45] as we 
found for Model I. This sharp fall can be understood in terms of the conflict between the 
light GP-even Higgs boson mass and muon g — 2 calculation. We see from the Aa^ — tan f3 
plane that there are plenty of solutions (green points) that satisfy all collider constraints 
and the requirement that 124 GeV < rrih < 127 GeV with tan f3 > 40. However, all of these 
solution provide negligible contribution to the muon g — 2 calculation. We conclude that 
for tan/3 > 40 we have heavy smuons or neutralinos in the spectrum in order to obtain 
the required Higgs boson mass. This observation can be related to the results presented in 
Ref. [44] , namely that the bottom squark and stau can provide negative contribution to the 
light GP even Higgs boson mass calculation. This contribution can become significant for 
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tan 0 


Figure 8: Plots in the Aa^ — tan/?, m/j — tan/?, m^o, and — tan/? planes. The color 
coding is the same as in Figure 1 except the Aa^ — tan/? plane whose color coding is the 
same as in Figure 5. 


tan/? >40. In order to compensate this and obtain correct Higgs boson mass, we need to 
increase the contribution from the stop loop. This, in turn, trough RGE evolution increases 
the gaugino and smuon masses. The mh — tan 0 plane conhrms this conclusion. Here we 
can see that it is possible to obtain the Higgs boson mass compatible with the experimental 
data (green points), with the muon g — 2 constraint presented in yellow points. 

Exploring the m^o — tan /? and — tan /? planes we found that the arc shape (green 
points) in these plots is caused by the Higgs boson mass requirement. Here we can see 
that the neutralino and smuon masses drastically increase for tan/? > 40, which leads to a 
correponding decrease in the supersymmetric contribution to the muon g — 2 calculation. 

Figure 9 represents the sparticle spectrum in the rriA — tan 0 and — rrig planes. 
The color coding is the same as in Figure 1. The — tan/? plane shows that can 
be as light as 1 TeV or so. This bound is also compatible with the muon g — 2 condition. 
Similarly the — rrig plane reveals relatively light squark and gluino masses. According 
to our results, muon g — 2 requires < 2.5 TeV, while rrig < 2 TeV. In contrast to the 
usual GMSB framework, the hrst two family squarks are lighter than the third family and 
quite possibly accessible at the LHG 13. 

Figure 10 displays our results in the Aa^ — plane. The color coding is the same as 
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Figure 9: Plots in the — tan/3 and — rrig planes. The color coding is the same as 
in Figure 1. 
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Figure 10: Plots in the Aa^ — and ttiq — planes. All points are consistent with 
REWSB. The color coding is the same as in Figure 5. In addition, the blue points represent 
the solutions with neutralino NLSP. The ttIq — rrig^ plane plots only the solutions with 
smuon NLSP. The color coding used in this plane is same as in Figure 1. 
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in Figure 5 and all points are consistent with REWSB. The blue points represent solutions 
with neutralino NLSP. As seen from the Aa^ — mg plane, the gravitino mass can he in 
a wide range from a few eV to 10 GeV, consistent with the muon g — 2 constraint. The 
blue region indicates solutions with neutralino NLSP which are not constrained by the 
prompt smuon decay. The residual green region is the set of solutions with smuon NLSP. 
The results for the gravitino mass are similar to those obtained in Model I that a gravitino 
mass lighter then 30 eV is only realized if the smuon is NLSP. However, in Model II, the 
solutions with neutralino NLSP , that are consistent with muon g — 2 anomaly yield a 
gravitino with mass of the order of GeV, which makes it difficult to have the latter as dark 
matter particle. 

Finally, we present two benchmark points in Table 3 to exemplify our Endings for Model 
II. All points are chosen to be consistent with the experimental constraints and muon g — 2. 
All masses and scales are given in GeV units except for the gravitino mass which is given 
in eV. Point 1 displays a solution with the largest tan (3 value compatible with muon g — 2. 
The NLSP happens to be the right-handed smuon for Point l.The staus are heavy, which 
push the right-handed smuon masses to low values through RGEs (in contrast to the left- 
handed smuons). Hence the large tan f3 region mostly yields right-handed smuon NLSP as 
displayed in the plots above. In addition. Point 1 represents a solution with the gravitino 
mass about 19 eV. Similarly, Point 2 represents a solution with gravitino mass of about 30 
eV. The spectrum is quite similar to Point 1, and one can compare how a slight change in 
the stau mass can afiect the right-handed smuon mass. The stau mass is slightly lighter in 
Point 2, which leads to a slightly heavier right-handed smuon with mass about 255 GeV 
in the spectrum. The benchmark points in Table 3 exemplify the impact of heavy staus in 
the low energy phenomenology. 


5 Conclusion 

We have explored the sparticle and Higgs phenomenology of a flavor symmetry based 
MSSM in the framework of gauge mediation supersymmetry breaking scenario. Explicit 
ultra-violet completion of models that generate a symmetry based MSSM (sMSSM) spec¬ 
trum [15] at low energies have been presented. These include models based on an additional 
U(l) flavor symmetry which acts only on the third generation fermions. The charge as¬ 
signments of these models (Type I & 2) are given in Table 1. Our model contains one 
additional parameter compared to the canonical GMSB scenario. Specifically we have a 
common SSB mass term for the first two family sfermions (which are lighter in mass com¬ 
pared to the third family sfermions.). This splitting between the family masses is generated 
by exploiting the U(l)symmetry which only acts on the third generation. This helps us 
to reconcile the muon g — 2 anomaly with the Higgs boson mass, and we can also accom¬ 
modate the desired relic abundance of dark matter in the form of gravitino. We find that 
the simultaneous explanation of these observables requires the gluino and first two family 
squark masses to be less than 2.5 TeV. The mass of the first two family sleptons is less 
than 600 GeV (or so) in Model I and less than 300 GeV for Model H. 


17 



Point 1 

Point 2 

A 

0.24 X 10® 

0.24 X 10® 

Afniess 

0.32 X 10® 

0.57 X 10® 

tan (3 

39 

38 

D 

11100 

10060 


13.7 

12.5 


4.9 

4.6 

Aa^ 

20.8 X 10-“ 

21.0 X 10-“ 

g 

4437 

4347 

At 

-334 

-350 

Ah 

-351 

-371 

A, 

-39 

-41 

rrih 

125 

124 

rriH 

1143 

1260 

rriA 

1136 

1252 

mH± 

1146 

1263 

m^o 

X-1,2 

382, 743 

341, 667 

m^o 

A3,4 

4364, 4364 

4277, 4277 

m-± 

Xl,2 

744, 4304 

667, 4221 

m-g 

2058 

1868 


2441, 2338 

2350, 2246 

^h.2 

10785, 11126 

9705, 10066 


2443, 2339 

2352, 2246 

mr 

01,2 

11085, 11281 

10025, 10232 


688 

686 

rriii^ 

33152 

30021 


770, 246 

763, 255 

TOfi 2 

32938, 33182 

29797, 30046 

rriQ (eV) 

19 

33 


Table 3: Benchmark points for Model II. All points are chosen to be consistent with the 
experimental constraints and muon g — 2. All masses and scales are given in GeV units 
except for the gravitino mass which is given in eV. Point 1 displays a solution with the 
largest tan/? value compatible with muon g — 2. The NLSP happens to be the right-handed 
smuon for Point 1. Point 2 depicts a solution with neutralino NLSP. 
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